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Abstract 

, We construct fermionic Lagrangians with anisotropic scaling z = 2, 

the natural counterpart of the usual z = 2 Lifshitz field theories for 
scalar fields. We analyze the issue of chiral symmetry, construct the 
' Noether axial currents and discuss the chiral anomaly giving explicit 

t— I , results for two-dimensional case. We also exploit the connection be- 

' tween detailed balance and the dynamics of Lifshitz theories to find 

. different z = 2 fermionic Lagrangians and construct their supersym- 

CNl ' metric extensions. 

• r-j ■ 1 Introduction 

X: 

Lifshitz field theories exhibit an anisotropic scaling between space and time 
that breaks Lorentz symmetry. Introduced in the context of condensed 
matter phenomena Q~|-[l], they are also actively investigated in connection 
with problems in gravity [5]-[7] and in particle physics (see [8] and reference 
therein) . 

The anisotropic gauge scaling for spatial coordinates x and time t takes 
the form 

x — > bx, 

t b z t (1) 



(N 



'Associated with CICBA 



1 



where b € R is the scaling factor and z G Z the degree of anisotropy which in 
condensed matter systems is identified as a dynamical critical exponent. The 
z = 1 case corresponds to the ordinary scaling law in the conformal group. 
In connection with this, it becomes natural to study, within the context of 
gauge/gravity correspondence, metrics with anisotropic scaling which could 
be dual to Lifshitz field theories in the same way as AdS metrics allow to 
establish duality leading to strongly coupled conformal field theories [5]. 

A prototype of a Lifshitz field theory is a (d+l)-dimensional scalar theory 
with dynamics governed by the action 



Here V 2 = didi is the spatial Laplacian. The scaling dimensions are [x 1 ] = 
— 1 and [t] = —z, with [0] = (d — z)/2 and [k] = z — 2. 

Models with action ([2]), particularly in the z = 2 case interest condensed 
matter physicists because they exhibit a line of fixed points parametrized by 
n and reproduce the phase diagrams of known materials [2]. Other z / 1 
bosonic field theories have been discussed, and also their gravity duals have 
been investigated using the holographic principle with a bulk geometry given 
(in Poincare coordinates) by the metric 



which for z = 1 becomes the usual AdS metric. 

Defining fermionic Lifshitz theories is an open issue. Already in refs.[3j|l], 
dealing with conformal quantum critical points, the problem is addressed for 
the case of z = 2 scaling by defining a ground state wave functional for the 
Hamiltonian of a 2 + 1 system in terms of the (1 + 1) Euclidean Dirac action, 
but the discussion in this work is restricted to a bosonized version of the 
theory. Self-interacting fermionic Hamiltonians leading to z = 2 dynamical 
critical exponent have been constructed in |10j in 1 + 1 dimensions following 
an approach different to that based in theories defined as in ([2]). More 
recently fermion theories with z = 3 have been studied in connection with 
its renormalizability, chiral anomaly and renormalization group flow |ll|-|17j. 

In this work we undertake the construction of z = 2 fermionic Lifshitz 
theory following a different way. Having in mind that the Dirac equation can 
be seen as the "square root" of the Klein-Gordon one, we firstly start from 
the bosonic action ([2]) and derive a 7-matrices algebra which naturally leads 
to a Lifshitz fermion model in d + 1 space-time dimensions. In doing so, we 
introduce the square root of the Laplacian as the source of the anisotropic 
scaling. We consider the z — 2 case which corresponds to (—V 2 ) 1 / 2 but the 
power z/2 could be equally considered giving rise to theories with arbitrary 
even z. We shall also discuss the issue of chiral symmetry for massless 




d d xdt [[dt^f - k 2 4> {-V 2 ) z 0) 
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fermions and consider the case in which they couple to an Abelian gauge 
field. We work out the corresponding classical Noether currents and we also 
compute the chiral anomaly in 1 + 1 space-time dimensions which, as it is 
well-known, is closely related to bosonization, one of the issues of relevance 
in condensed matter applications of Lifshitz theories [3] . 

Secondly, we present an approach based on the connection between the 
principle of detailed balance and the dynamics of Lifshitz theories stressed in 
refs. [5]-[7]- The main ingredient in this approach is stochastic field theory 
quantization provides, a useful tool in the construction Lifshitz field theory 
actions. Indeed, starting from the relativistic action for scalar fields in d 
Euclidean dimensions stochastic quantization leads to a d+ 1 effective action 
with the fictitious time derivatives appearing quadratically and those of the d 
remaining derivatives quartically. Following this approach, we shall analyze 
the stochastic quantization of fermionic theories and construct various Dirac- 
Lifshitz actions with different anisotropy scaling z including z = 2. 

The paper is organized as follows. We establish in Section 2 a connection 
between z = 2 bosonic and fermionic theories leading to a consistent z = 2 
Dirac-Lifshitz action in d+ 1 space-time dimensions. In section 3 we discuss 
the issue of chiral symmetry finding the classically conserved axial currents 
an calculating the quantum chiral anomaly in the particular case of a 1+1 
dimensional space. Working within the stochastic quantization approach, we 
discuss in Section 4 the construction of alternative z = 2 fermionic theories 
in d + 1 dimensions as well as their supersymmetric extension. Finally we 
present in Section 5 a discussion of our results. 

2 Connecting bosonic and fermionic z = 2 Lifshitz 
theories 

Consider a Lifshitz bosonic action of the form ([2]) in d+1 Euclidean spacetime 
with a z = 2 scaling 



S L 
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...,d (5) 



It will be useful for our purposes to write action ([5]) in the form 





with G av defined as 
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The integral on the r.h.s. of eq.Q can be regarded as a degenerate scalar 
product on the space of vector fields on M, 

( ,) M :X(M)®X(M) — >-R 

such that 

(V(x), W{x)) M = j d d x (V W - K 2 ViV 2 Wi) (8) 

The associated group of isometries can be obtained from the invariance re- 
lation 

(A(V),A(W)) M = (V,W) M 

which explicitly reads 

A 0lM A 0u - K 2 A ilM A iu V 2 = (V 2 ) (9) 
Taking /z, v = 0, we get 

A i0 Aio = , Ag = 1 

which implies 

A i0 = , A 00 = ±1 

respectively. Alternatively by taking \i = i and v= j in equation ([9]), we 
determine the relations defining the remaining components, namely, 

Aoi = A^Akj = 5{j 

implying that Aj,- are the entries of an orthogonal matrix. In summary, the 
final form of A is 

A 00 = 1 , A 0i = A i0 = , A ki A kj = ff y (10) 

meaning that that isometries of G^y (V 2 ) reduce to rotations in the spatial 
sector and time reversal map. 

One can extend this framework to the case of fermionic theories in a 
natural way by defining "local" gamma matrices T satisfying 

{r^, r^} = 2G I1U 

An appropriate choice for these gamma matrices is 

Tt = K (-V 2 ) 1 / 2 7i , r = 7 o , i = l,2,...d 
where ( 7a ) = ( 7 o,7i) satisfy the ordinary Euclidean-space Clifford algebra 
{7a, 7&} = 2 $ab, a = 0,l,...d 



4 



and (—V 2 ) 1 / 2 is the nonlocal operator which can be defined as in |18| (see 
also references |19|-|22|). 

{-V^f{x) = P.V.I HX) ~ ^ d d y (11) 

jR d \x — y\ 

Defined in this way the square root of the Laplacian is an elliptic pseudod- 
ifferential operator of order 1. However, a better insight of the non local 
operator (—V 2 ) 1 / 2 is attained in the context of the harmonic extensions of 
a compactly supported function on R n to M. n x (0, oo), as in ref. |18| . where 
it appears as a map from the Dirichlet-type data problem to the Neumann- 
type data one. Moreover, it can be shown that it is self-adjoint and positive- 
definite (see |23|-|24| and references therein). In Appendix we give a brief 
description of this approach and some properties of this operator. From now 
on, we shall use a more compact notation for the square root of the Laplacian 
naming it by T 

T : = (-V 2 ) 1/2 (12) 

For further developments, it is useful to work out the action of the operator 
T on imaginary exponentials. From the identity 

— V 2 exp(z/cx) = | k\ 2 ex.p(ikx) (13) 

and the natural definition of the square root Laplacian it can be established 
that 

T exp(ikx) = \k\ exp(ifex) (14) 

(see the Appendix for a proof of this result). Let us point that, strictly 
speaking, since the operator symbol \k\ is not smooth at the origin, the 
problem should be handled inserting an appropriate cutoff function. 
We can now propose a vierbein-like connection between V and 7 

(V) = (o ® ti ) S(x - y) , (£/)=( _I TJ ( 15 ) 

lowering and raising indices with G^ u and 5 ao so that 

(E^ a E ai/ ) = (G^y) , (E a u E ub ) = (6 a b) 

with products including an integration over space-time. 

This set of matrices {r^} naturally induces a definition of a Dirac-Lifshitz 
equation in the form 

(iT^ - m)V = (16) 

Here m has dimensions of mass squared. In terms of the ordinary 7^ matrices 
eq, ([T6|) reeds 

{ijodo + inT(jidi) - m) ^ = (17) 
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This is nothing but the Euler-Lagrange equation for action 

Sdl = j d d+1 xi> (i l0 d + mTjA - m) V (18) 

These Dirac-Lifshitz fields shares the anisotropic gauge scaling ([1]) for z = 2. 
From this expressions we define the Dirac-Lifshitz operator 

7^ = ij d + inT^idi -m (19) 

It is interesting to note that, as in the ordinary Dirac equation case, 
one can proceed to "square" the Dirac Lifshitz equation looking for a result- 
ing Klein-Gordon-like equation. Indeed, multiplying (|16p equation by the 
operator (iT^d^ + m) and using the identity (see Appendix) 

Tdi = 3iT (20) 

we get the usual z = 2 Lifshitz equation for scalar fields, subject of many 
recent investigations 

(d 2 - K 2 (V 2 ) 2 -m 2 )^ = (21) 

This is precisely the equation of motion arising from the z = 2 bosonic 
action ([!]) and in this respect Lagrangian (fT8|) seems to be a natural candidate 
to study Lifshitz fermionic models in connection with the problems of the 
conformal quantum critical points [3|- 

3 Chiral symmetry and the z = 2 Dirac-Lifshitz La- 
grangian 

The issue of chiral symmetry and the quantum anomaly of z = 3 Lifshitz 
theories with massless fermions coupled to a gauge field in 3 + 1 space- 
time dimensions has been thoroughly studied in [12], |15|-|16|. Using the 
path-integral framework the resulting anomaly A was shown to be identical 
to the usual relativistic calculation and this is likely to be related to its 
topological character in the sense that A is metric indepedent and its integral 
-the topological charge- does not depend on local details but only on global 
properties. The calculation for the z = 2 Lifshitz defined in the precedent 
section seems more involved since it implies handling the square root of the 
Dirac operator both in the action and as the natural regulator in a sensible 
definition of the path-integral measure. We shall now discuss this issue. 

3.1 The free massless fermion theory 

In the m = case the free fermion Lifshitz action (|18p is invariant under 
chiral rotations 

ip -> exp(7 5 e)V> , i> ->■ V5exp(7 5 e) (22) 
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There should then be, at the classical level, a continuity equation for the 
chiral current and an associated conserved chiral charge. To obtain their 
explicit form one can proceed as follows. One promotes the global change to 
a local (infinitesimal for simplicity) one 

tP ->. if/ = exp(e(x)7 5 )V> ~ (1 + e(x)^ 5 )ip 

-> ? = $exp(e(x) l5 )^ij(l + e(x)^) (23) 

under which the action should necessarily change, 

5S DL = J d d+1 x$(l + 75 e(x)) (i 7 °5 + iKTfdi - m) (1 + -y 5 e(x))^ - S DL 
After some work one finds 

5S DL = j d d+1 xi^ 1 j 5 ^d Q e + j d^xin^T^^die 

= -i f d d+1 x^ (7°75^5 + kTY 75 ) (24) 



For constant e this term should be absent and then we get a continuity 
equation and the explicit form of the chiral current 

« = o 

it = i>i°ihi> , 3% = KtPiVip (25) 

As it was to be expected, the chiral charge coincides with the usual (Dirac 
theory) one, 

Q s = f cf^ 7 °75^ = / d d x^^ (26) 



There is of course a £7(1) global invariance with associated current 

i = $yV, f = K^T^ (27) 

so that again, the conserved £7(1) charge coincides with the usual one, 

Q = [ d d x^°tp = [ d d x^ij (28) 



Coupling fermions to £7(1) gauge fields 

In order to couple massless z = 2 Lifshitz fermions to a £7(1) gauge field 
background A^, we introduce the extension of the Lifshitz-Dirac operator 

m 

V = ^V^A] = 7oD [A] + KrfxDx [A] T [A] (29) 
where we have defined 

T[A] = (D j [A)D j [A}) 1 / 2 (30) 



7 



with covariant derivatives given by 



D = ido + eA , D { = idi + eAi (31) 

Some properties of operator T[A] are discussed in the Appendix. In par- 
ticular, T[A] is an Hermitian operator so that T> defined as in (|29|) is also 
Hermitian. 

A z = 2 Dirac-Lifshitz action can be compactly written in terms of the 
Dirac-Lifshitz operator as 

Sdl[A] = J d d+1 x4,^V^[A]^ (32) 

As in the free case, this action remains invariant under global chiral transfor- 
mations (as well as under local U(l) gauge transformation). So, there should 
be a classically conserved axial current. Concerning the quantum level, one 
should expect an anomaly which can shall calculate within the path-integral 
framework following the well-honored Fujikawa's approach |25|-|26], We start 
from the partition function written in the form 



/ 



D^Di/}exp(-S[A]) (33) 



and consider the infinitesimal local path-integral change of variables given 
by eq. (|23p . The action in the exponential now gets, as in the case of the 
free fermion classical action, a contribution proportional to d^e; and con- 
cerning the path-integral measure, there is a Fujikawa Jacobian Js[e] . So, 
the partition function changes to 

Z = J Dj/Di//J B [e] exp (- S[A] - 5S[A; d^eU (34) 

As in the free fermion case, we shall show below that 5S can be accommo- 
dated as 

6S[A;d ll e] = J d d x(dA*M (35) 
so that Z can be written as 

Z = J Dj/DtfJ 5 [e] exp (- S[A] - e(x)(d^)) (36) 

Since the generating functional cannot depend on the parameter e(x) intro- 
duced by the change of variables one has 

(37) 

which leads to the anomaly equation 

M)=A (38) 
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where 

oe[x) 

Let us now prove eg. ([35]) . The action S[A] can be written as 

S[A] = J d d+1 xij ( 7o A) + K-f (Di[A]T[A])) $ (40) 

After the change of variables (|23p one has 

SS[A] = S [A] + Si[A]-S[A] (41) 

with 

S [A] = yd d+1 ^(l + 75e(x))7 D (l + 75e(^)^ 

SM] = k /d d+1 ^(l+ 75 e(x))7 i Ar[A](l + 75e(x))^ (42) 



The term So contains the same variation arising in an ordinary Dirac theory 
and can then be written in the form 



S = - I d d+1 xe(x)di (V> 7 °75^) + / d d x&j°Do4> (43) 
In order to handle 5S\ in a similar way, we first integrate by parts 
SS 2 [A] = k f d d+1 xD*^{l + 7 5 e(x))7 i r^] (l + 75 e(x))V 

— k j d^xijjfDiTiA}!/} (44) 

finally getting 



SS 2 = -in J d d+1 x($e)V'7VT[A]V> = in J d d+1 xe(x)di (^VTLA]^) 

(45) 

Putting all this together we have for the complete variation of the action 

SS = j d d+1 xe(x)d^ (46) 

with 

i 5 ° = ^ 7 V, ji = ^YT[A}4> (47) 

which differs from the free fermionic case because j§ contains the operator 
T[A] = (AA) 1/2 instead of T = (-V 2 ) 1 / 2 and hence depends on A { . 

We have now to compute the Fujikawa Jacobian associated to the chiral 
change of variables. For simplicity we shall consider the 1 + 1 dimensional 
case in which Lagrangian in action (J32]) takes the simple form 

L = ^°(D )ip + K^(j 1 D 1 r[A)) V> (48) 
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Our choice of gamma matrices is 

7o= (io)' 7l = (-*o)' t* = **™ = ( o -°i) (49) 

Following Fujikawa's approach, we consider a change of variables as in 
(|23p and calculate the associated Jacobian J 5 [Je(x)] which gives 

J 5 [Se(x)\ = exp ^2 J d 2 x ^ ip\ l {x)^ n {x)8e(x)^ (50) 

From the fact that the generating functional cannot depend on the parameter 
e(x) introduced by the change of variables 

62 (51) 



Se(x) 

one then gets 



5e(x 

with A given by the ill-defined formula 



a ^ = S ^ =A (52) 



A = 2'Vvl(xHMx) = lim2tr 75 «< 2 l(i-i;) (53) 

n 

which thus requires an appropriate regularization 



Area = Hm 2tr7 5 5 (2) (x - y) 

y^tx 



(54) 



leg 



A practical regularization approach is the so call heat-kernel method 
which consist in introducing an appropriate regulating operator R 2 and a 
cut-off mass M in the form 

Areg = M 2tr ( 75 exp(-i? 2 /A[M] 2 )5 {2) (x - y)) (55) 

where A[M] = M q with q chosen so that [A] = [R] = q. A finite result for 
the anomaly is obtained by interchanging trace and limit 

A= lim lim2trf7 5 exp(-i? 2 / A2 )^ (2) (x-y) N ) (56) 
M^ooy— >x \ / 

Now, the regulator should be chosen so as to preserve the basic invariances 
of the theory, in the present case gauge invariance and also the invariance 
under the anisotropic scaling ([TJ . This makes natural to use as regulator the 
operator appearing in the Lagrangian so that R? = T>T> and then one has to 
choose A[M] = M 4 . 
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In two space-time dimensions the Dirac-Lifshitz operator reduces to 

V = i lQ D + k^Di[A]T[A] (57) 

and T[A] is simply 

T[A\ = {D l [A\D 1 [A]) 1/2 (58) 
In the Appendix we show that T[A] can be also written in the form 

T[A} = W~ 1 [A 1 }TW[A l } (59) 

where 

W[Ax] = exp f-ie J A 1 (y) dyj (60) 

At this point, it is important to note that one can gauge away one component 
of the gauge potential just by a gauge transformation. In particular, we shall 
consider WL4i] as a ?7(l)-gauge group element W [Ax] so that 

(An, Ai) — > (W [Ax] AqW- 1 [Ax] , 0) 

which is formally implemented as 

Ap—tAp + W [Ax] (ify) W' 1 [Ax] (61) 

This gauge will greatly simplify calculations leading to a result that can be 
trivially extended to a general gauge. 

The regulator T> 2 can be now written in the form 

V 2 = ID D -iK 75 [D^W- 1 [Ax] (idjTW^+tflW- 1 [Ax] (-$) 2 W[A 1 ] 

(62) 

Under the gauge transformation (|61jl D„ changes according to 



D„ — > W [A,] D^W- 1 [Ax] 
so that the regularizing operator becomes 

W [Ax] V 2 W~ X [Ax] = IW[A 1 ]D D W~ 1 [A 1 ] + k 2 I (-d 2 ) 2 

-m 75 [W [Ax] DqW' 1 [Ax] , (idx)T] (63) 

with this and the 5^ (x — y) function written in term of the plane wave basis 
eq. ([56"l) becomes 



A * Um Um f £ p2 tr Ue-ip^eHw^W-HA^/M^ip^ 

9 (2vr) 2 M^oov^xJ F \' a 

(64) 

Taking now the limit y — )■ x, the operator content of A reg becomes 
e -ip^ W [Ax] V 2 W~ l [Ax] e ip ^ = J$ (A,p) - zre 75 (A,p) + p 2 + K 2 Ip\ 
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Here we have used that 

and introduced <3? and defined as 
<f>(A,p) = W 



W 



A, + ^ 

e . 



A 1 + *: 

e 



Q(A,p) 



W 



A l + P -± 

e 



D 2 - 2 P0 (id )) W- 1 
(id ) W- 1 



A l + P -± 



A\ + ^ 

e 



(idi)Te ipiXl 



A 



Making the substitution po = M 2 k§ and p\ = Mk\, the anomaly becomes 
1 



leg 



lim J d 2 p 2M 3 (exp - (jfejj + K 2 kf)) 



xtr ( 75 exp(-(/$ (A,M 2 A;o,MA;i) - iK~( 5 G (A, M 2 k , Mh)) /M 4 )) 

(65) 

By M power counting in the argument of the exponential, one concludes that 
only the linear term in G [A, M 2 ko, MkA survives after taking M — > oo 
and the 7- matrix trace. In order to carry out the explicit calculation of 
this contribution, recall that before differentiating with respect to e{x) the 
regulated Jacobian includes an integral over 2-dimensional space and hence 
one can use the integration by parts properties of T (see the Appendix) and 
then differentiate, thus obtaining 

Q [A, M 2 k , Mki) =ieM\ki\F 01 

The anomaly can then be finally written in a gauge invariant way as 



A 



reg 



4vr £ ^ 



(66) 



which coincides with the result for the ordinary relativistic two dimensional 



-— F 



(67) 



As observed in refs. |12| and |15j for z = 3 models in 3 + 1 dimensions (and 
then generalized to arbitrary values of z in |16j). the coincidence between 
the anomaly of relativistic z = 1 and Lifshitz z / 1 theories is not surprising 
in retrospect since the anomaly A is related to a topological density that 
coincides with a total divergence, so that eq. ()67p is the only possible result 
for the anomaly, including the universal coefficient, independently of the 
value of z. To confirm this en the 1 + 1 dimensional case, we briefly describe 
the 1 + 1 anomaly calculation for the z = 3 model, with the action in the 
form usually taken in the literature, 



S 



(68) 
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Here the Dirac-Lifshitz operator is given by 

2,(3) =7 M2,(3) = 7 o Do + ly {Di{ _ DkDk) + (- Dk D k )Di) (69) 

with the covariant derivatives defined as in (|3ip . Again, at the classical the 
there is a conserved Noether charge associated to a global chiral transforma- 
tion. At the quantum level, one should take into account the change in the 
fermionic measure and one ends with the anomaly equation (|52p 

A= lim lim2tr f7 5 exp(-P (3) P (3) /M 6 )5 (2) (x-y) N ) (70) 

As before, the regulator has been chosen so as to preserve the basic 
invariances of the theory (i.e. gauge invariance and also the anisotropic 
scaling invariance) and hence A[M] in (|55p should be A[M] = M 6 . 

In 1 + 1 dimensions the Dirac-Lifshitz operator reduces to 

V® = i 7 oA) - KjiDf 

and X?( 3 ) 2 is given by 

V^ 2 = ID 2 - ivys [D ,Df] + k 2 ID\ (71) 



We shall not repeat the details of the calculation that follow closely those 
that leading to (|6T[) in the z = 2 case. We just quote the only non-vanishing 
M — > oo term in A which instead of integral (|67p is now given by 




which again coincides with the ordinary result for the two-dimensional anomaly, 

4 Lifshitz fermions and stochastic quantization 

It has been recenlty observed that field theory models with the anisotropic 
scaling ([T]) can be thought of as stemming from a stochastic quantization 
process [6]- [7], |28]-|29|. In this way, a d + 1 dimensional bosonic Lifshitz 
theory like that defined by action ([3]), with z = 2 can be thought off as 
descending from a free bosonic action in d Euclidean dimensions, 

= I d d xd^d^ , i = l,2...,d (72) 
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The dynamics in extra time variable r is ruled by a Langevin equation, 
namely a first order diferential equation with a dissipative (drift) term and 
driven by a source of stochastic noise. As far as the quantization of a d- 
dimensional theory is concerned (the asymptotic equilibrium state), there 
exist some arbitrariness in the election of this stochastic process. One may 
take some profit from this freedom to build up theories in d + 1 dimensions 
with different behaviours under "time" and space scaling. 

In the case of the bosonic field (j3J), to get a z = 2 Lifshitz theory one 
introduce an extra dependence of the fields on this extra time variable, 
<p(x) — >■ (/)(x,t), and compose the Langevin equation with 6S^ /6(p playing 
the role of a drift force, and a white noise source as the stochastic force |30j. 
The case of fermionic fields can be treated analogously |31|-|32j. with ip and 
ip taken as independent fields, each one with its own Langevin equation 



dr 6ip 
dtj>(x,r) 5S[i/),il>] 



+ fj(x,r) (73) 



dr 5ijj 
where rj, fj are anticommuting fermionic gaussian noises satisfying 

(Va(x,T)) = (fj a (x,T))=0 (74) 

{r) a (xi,T 1 )fjp(x2,T2)) = 26 (h -t 2 )6 ( - d \xi - x 2 )6 a( 3 (75) 

with a, j3 the spinor indices. Such normalization arises from the following 
definition of the (Fokker-Planck) generating functional 

Z[J,J] = J DfjDrj exp (— - j dx j dr'fjT]] exp (J dx j dr 1 (fjJ + Jri)\ 

= exp ^2 J dx J dT'JJ^j (76) 

where spinor indexes have been omitted and J, J are fermionic external 
sources. 

Now, as already noted in [31j-[32j, when S is taken as the Dirac action 
one should in general include an appropriate kernel K(x, y) to ensure that 
the relaxation process is such that the system tends to the equilibrium when 
r — > oo. Different kernels produce stochastic process which can differ in 
the behaviour under scaling. In any case, Langevin equations (j73|) are then 
replaced by 

mX ' T) ~ f'dyK(x,y) 6 -^4 +V (x,r) 



dr J ' 6ip(y,r) 

f dy 6 41M^) + v( X ,r) (77) 



dr J 6\fj(y,T 



14 



and eq. ([75]) by 

{Va(xi,T 1 )fj/3(x 2 ,T 2 )) = 28{ti -t 2 )K{xi - x 2 )5 a p (78) 
The usual choice for K 

K(x,y) = (i d x + m)5(x-y) (79) 
in the case S is the Dirac action for a massive fermion, 

S= [ d d xip(x)(i @-m)ip(x) (80) 



gives rise to the Langevin equations 



^ (X ' T) = (V 2 -m 2 )iP + ri(x,T) 



dr 

(V 2 -m 2 )i; + fi(x,T) (81) 



dip(x,r) . 2 



dr 

where a z = 2 anisotropic scaling turns evident. 

Let us now come back to the general fermionic stochastic process defined 
by the Langevin equations ([78]) . In the presence of the kernel K the partition 
function associated to (|76|) takes the form 

Z = jDfjDriexp(-^ j ' d d xd d ydTfj(x,T)K-\x,y)r)(y,T)\ (82) 

We now proceed to change variables from fj, Tj to tp, tp 

Z= f D^DVdet" 1 ^det" 1 exp (-- [ d d xd d ydTf](x, t)K~ 1 (x, y)ri(y, r) 
J dip dip V 2 J 

(83) 



Using Langevin equations (|8ip written in the form 

K i -- = = + K L rj 

dr Sip ' 



we end with 



x exp| 



dr 5ip5ip J \ dr dipdip 

l - fd'xdrl^K- 1 &S[M] \( ty{*,T) ]K SS[^] 
2 J \dr dip )\ dr Sip 

(85) 



15 



where we have used the identity 



' K-]d T ^ = K- a ^d T (86) 



Determinants in (|85p can defined as the (regularized) product of eigenvalues 
of the equation 

((*-r^-^)^? (87) 

and a similar one for the other factor. These determinants can then be 
represented as quadratic path-integrals over two complex commuting fields 
4>f (i = 1,2) carrying a spinorial index which can be identified with the 
spinorial ghosts for ghosts |33| . With this, the partition function Z can be 
written as 

Z = J DtfD^DfaDfaDfaDfaexp^- J ' dxdrL eff (88) 
where 

1 f d^(x,r) x SS[j^\\ f d^(x,r) „6S$,4] 
L eff = o 7^ K EZ ^ + K 



2 \ dr 5ip J \ dr Sip 



dr 5ip5ijj J \ dr 5ip5ip 



92 

(89) 



or 



_ 1 %,T) ^(1,7) l ^feVj y ggfe^ 

e// " 2 dr Or 2 5^ 6$ 



(90) 

We can see now how the choice of the stochastic process, defined by a kernel 
K, gives rise to different scaling behaviours. Indeed, starting with a d - 
dimensional free massive fermion Dirac action and using a trivial kernel 

K(x,y) = 6(x-y) (91) 

(this being consistent for m ^ the fermionic part of the effective fermion 
Lagrangian (|90jl will have second derivatives both in time and space variables 

tM = f*£ T)8 *£ T) +l(VW))» (92) 
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which corresponds to a Lifshitz theory wit z = 1, i.e a theory with isotropic 
scaling. If one instead chooses as d dimensional action a Klein-Gordon 
fermionic action 

S = [ d d x$(x)(V 2 + m 2 )^{x) (93) 



one ends with a z = 2 Lifshitz theory with Lagrangian 

J 2 OT OT 

We shall end this section by making use of another remarkable attribute 
of the stochastic quantization scheme, namely the inherent supersymmetry 
of the partition function, in order to supersymmetrize the Lifshitz fermionic 
models constructed above |34] . Alternatively, and following the ideas intro- 
duced in the section 2, we may consider the d dimensional action for fermions 

Slf = f d d xijj (^(-V 2 ) 1 / 2 ^ - m) ij> 

which gives rise through the stochastic quantization approach to the fermionic 
Lifshitz Lagrangian 

= ft T> ^ - ¥ {-^-^iM-vy>hA + ^) * 

= l«^MEl[) + I;, (kW - rrf) i, (95) 

which also exhibits a z = 2 anisotropic behaviour under scaling. 

Starting from partition function Z as given in (|88p with Lagrangian (f90|) . 
we introduce fermionic auxiliary fields F,F so that Z becomes |35j 



DipDipDfaDfoDfaDfoDFDF exp ^- ^ d d xdTL eff ^,ip,(pi,(f>2,F,F] 

(96) 

We now define fermion superfields in the form 

^ = rp + Ofa + fce + mF 

* = $ + efo + (j)i9 + i90F (98) 
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and supercovariant derivatives D and D acting on superfields according to 

d d - d 

d 

D 2 = D 2 = , {D,D} = — — = —H (100) 
Supersymmetry generators Q and Q take the form 

°-| +ff 5F (101) 

so that 

{Q,Q} = J^ = # (102) 

A supersymmetry transformation with parameters (e, e) on superfields reads 

Qm = eQ + Qe (103) 

In order to write the supper symmetry transformations in component 
fields we consider the action of Q and Q on superfields 

Q^> = -^Lu + Ofa + 4> 2 + i69F) = fa + i6F 
39 

= (Jg+ O^A (V + Hx + hfi + i96F) 

= 4> 2 - of + 9ip + ee$2 (104) 

so that 



and then 



Qm = e (01 + iOF) + ( <p 2 - iOF + 0if> + 99<p 2 ) e (105) 

(106) 



<50 = + <50i = e(ip — iF) 

54>2 = itF 5F = —i&e 

Let us consider the action 



Si = J d9d9D^Dt> (107) 
that written in term of component fields is given by 

dOdO— (tjj + 94)2 + 4>\9 + i99F) 

= J d9d9 (fa - i9F) (fa + i9F - 9ip + 09 fa 



h<Pi+ FF + iFip (108) 
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Define also 



So 



J d6d6D^Dty (109) 



or 



dOdO^ (ip + 0<j>x + foO + i66F) 
39 dr 



J d9d9 (4> 2 - iOF) (<fc + iOF -6^ + 99 
0202 + FF + iFip (110) 



so that 



St + S 2 = J d6d6 (D^D^ + DVDV) = 0i0i + 2 2 + 2FF + i{F^j + ^F) 

(111) 

For the case K = K 1 = I the terms in the right hand side exhaust all 
terms in (|97|) containing temporal derivatives. For nontrivial A one has just 
to make the following change of variables: 

FA" 1 ->• F , F -»• F 

TpK" 1 — > -0 , if) ^ Ip 
4>lK~ l -> 01 , 01 -> 01 

02/f" 1 -> 2 , 02 ^ 02 (H2) 

and the calculation reduces to the A' = A'" 1 one. 

As for the terms in (f97|) depending on action 5, one has just to add to 
Lagrangian ([Hip the following term 

S S = j d9d6Sfi>, *] (113) 

Expanding the action in powers of the superfield the only terms that con- 
tribute to the 9, 9 integral are 

S S _ / *»*[*. •] - F| + | F + + (114) 

so that the fermionic sector of superaction 

S = J d9d9d d xdr (D^D^ + DVDV + Sf[V, *]) (115) 
coincides with the action defined from Lagrangian (|90p . 
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5 Discussion 



The source of the anisotropic scaling in the Lagrangian formulation of Lif- 
shitz theories is the choice of spatial and time derivatives of different orders. 
This is generally achieved by replacing first order spatial derivatives in rela- 
tivistic Lagrangians by the Laplacian operator (or integer powers of it) while 
leaving unchanged the term with time-derivatives. Concerning bosonic the- 
ories, whose standard relativistic Lagrangian is quadratic in spatial deriva- 
tives, such replacement gives rise to a z = 2 anisotropy (or z > 2 for higher 
Laplacian powers [8]). For fermionic theories, application of the same strat- 
egy to the Dirac Lagrangian leads to an anisotropy degree z = 3 (or z > 3 
|16| ) . Now, Lifshitz theories have recently received much attention in connec- 
tion with condensed matter systems for which the z = 2 case is of particular 
interest. In particular, in [3], [I] the issue of a z = 2 fermionic systems in 
2 + 1 dimension has been related to a bosonic z = 2 Hamiltonian, invoking 
a bosonization approach. 

We have proposed in this work another route to study Lifshitz fermion 
actions with z = 2 in d+ 1 space-time dimensions, by introducing the square 
root of the Laplacian, and its natural extension when fermions are coupled 
to gauge fields. Roughly, it consists in replacing spatial derivatives d{ by the 
operator (— V 2 ) 1//2 <9j, which naturally leads both to bosonic and fermionic 
z = 2 Lifshitz theories. We have stressed (and briefly explained in the 
Appendix) that although (—V 2 ) 1 / 2 is a non-local operator, it can be handled 
as a local one in the context of the harmonic extension to an additional 
dimension. Exploiting this fact, one can prove useful identities that allow to 
discuss relevant physical properties of the systems with dynamics governed 
by Lagrangians containing such operator. 

A z = 2 Lifshitz action for massive free fermions was presented in eq.(|18|) 
and for massless fermions coupled to a U(l) gauge field in eq. (|32|). It should 
be pointed that "squaring" the resulting Dirac-Lifshitz equation (|17|) one gets 
the prototypical equation for Lifshitz bosons (eq.([2J) with z = 2) in the same 
way as the Klein-Gordon equation can be obtained by "squaring" the Dirac 
one. 

The z = 2 Dirac-Lifshitz action that we have proposed for massless 
fermions is invariant under chiral transformations and we have constructed 
the associated Nother current which leads to classically conserved charges. 
The j'5 component coincides with the usual chiral charge density arising in 
the relativistic case but the spatial components j$ differ. They have how- 
ever a very simple form given by eq. (|27p for free fermions and eq. (|47p for 
fermions coupled to an Abelian gauge field. This should be contrasted with 
the z = 3 case in which jf turns to be much more complicated |15| . In 
any case the coupling of gauge fields to Lifshitz fermions makes the resulting 
chiral current depend on the gauge field even in the U(l) case, in contrast 
with what happens with the relativistic theory. It should be noticed that 
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the same phenomenon takes place for the case of z = 3 theories constructed 
from local operators I16| . 

As it was to be expected, there is an anomaly at the quantum level which 
we have calculated in the two-dimensional case, showing that it takes the 
same form as in the ordinary Dirac theory. The fact that the anomaly is the 
same for relativistic and Lifshitz fermions is in agreement with the results in 
refs. |13]-|15tfTE] for the z = 3 case in 3 + 1 dimensions, and should be related 
to the connection between the obstruction to the chiral current conservation 
law and the topological charge of the gauge field configuration. 

In view of the connection between the anomaly in two dimensional space- 
time and bosonization (see |27] and references therein) one could expect to 
find the rules to connect the z = 2 fermionic theory we propose with a z = 2 
two-dimensional bosonic theory, possibly with action ([2]). We hope to come 
back to this problem in a future work. 

We discussed in section 4 fermionic Lifshitz theories with z = 2 within the 
stochastic quantization framework, a natural one recalling that in Parisi and 
Wu approach [30] one passes from a d dimensional Euclidean action to a d+ 1 
effective actions which in general exhibit anisotropic scaling. Since one can 
introduce in the Langevin equation different kernels acting on the drift force, 
one can find different effective d+1 actions which nevertheless correspond to 
the same d- dimensional quantum theory. Exploiting this fact we were able 
to construct several z = 2 fermionic Lifshitz theories. Moreover, through 
the well-known connection between the Fokker-Plank partition function and 
supersymmetry, we were able to find the supersymmetric extensions of the 
resulting fermionic z = 2 theories. 

Acknowledgements: We would like to thank Horacio Falomir, Eduardo Frad- 
kin, Gustavo Lozano, Mariel Santangelo and Luis Silvestre for helpful com- 
ments and suggestions. F.A.S is associated to CICBA and partially sup- 
ported by CONICET, ANPCyT and CICBA grants. H.M. thanks to CON- 
ICET for finantial support. 

6 Appendix: Square root operators 

In order to handle the operator (—V 2 ) 1 / 2 on R™ we follow ref. |18| where it is 
defined in terms of the harmonic extension problem to the upper half space 
R™ x (0, oo) of functions on M n . That is, given a function / : W 1 — >M, one 
looks for an harmonic function u : W 1 x (0, oo) — >M such that its restriction 
to M. n coincides with /. Thus, the problem reduces to solve the Dirichlet 
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problem 

u(x,0) = f{x) 

(116) 

V 2 n+1 u(x,y) = 

where V^ +1 is the Laplacian operator in M n+1 , x S M. n and t/£l. This is a 
well studied problem: for a smooth function on Co°(K n ) , there is a unique 
harmonic extension u £ C°°(lR n X (0, oo)). 

One then introduces an operator T defined on functions / : W l — >M. 
having an harmonic extension u : R n x (0, oo) — >M, such that 



(Tf)(x) = - 9u ^ y) 



Oy 



(117) 



Thus, since (Tf) (x) has also an harmonic extension to ~R n x (0, oo)), namely 
u y (x,y) provided u(x,y) is the harmonic extension of /, the sucessive ap- 
plications of T gives 



Oy 

that is equivalent to write 

T 2 f(x) = (-V 2 n )f(x) 

Then one can identify 

7 2 \ 1/2 



= - V 2 n u(x,y)\ =-V 2 J (x) 



T=(-V z n ) (US) 



Let us consider the case n = 1, so / : R — >M. and we look for harmonic 
extensions to K x (0, oo). For instance, taking / (x) = coskx, it is easy to 
see that u (x, y) = e~^ y cos kx, y 6 (0, oo), is a bounded harmonic extension 
of fix), 

+ q~2~J U ( X ^) = -k 2 e~\ k \ y coskx + \k\ 2 e~ |fc|s/ cos kx = 

Hence, we may evaluate the action of the operator (— Vf) 1//2 = (— d 2 ) 1 ^ 2 on 
/ 

(— d 2 ) 1 ^ 2 cos kx = T (cos kx) = — u y (x, y)\ y=0 = \k\ cos kx 

The same arguments hold if one takes / = sin kx since the function v (x, y) = 
e> k ' y sin kx is a bounded harmonic extension onlx (0, oo) of sinfcx 

( ^2 + J ^(^'2/) = -^ 2 e~ |fc|y sinfcx+ |A;| 2 e" |fc|?/ sin kx = 

and 

1 /2 

(— (9^) sin kx = T (sin for) = — v y (x, y)\ = \k\ sin fca; 
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We can then write the action of (— df} 1 ^ 2 on e tkx as 

(-df) 1/2 e ikx = (-d 2 ) 1/2 coskx + i(-dj) 1/2 smkx 

or 

(-d 2 ) 1/2 e ikx = \k\e ikx (119) 

this proving eq.(16). 

Remark 1: Suppose that the function / (x) admits an harmonic extension 
u(x,y), so (Tf)(x) = —u y (x,0). Then, functions dif (x), i = l,...,n, have 
harmonic extensions diu(x, y) and 

(T ($/))(*) = -d y (d iU (x,y))\ y=0 

= -diu y (x,y)\ y=0 

= di(Tf)(x) 

provided that u £ C 2 (R x (0, oo))). From this, we conclude that 

(-v*) 1/2 a = fl! [ (-v*) 1/a (120) 

Remark 2: We now address the selfadjointness of the operator T, following 
reference [24j. Let u(x,y) and v(x,y) harmonic extensions to C = M. n x 
(0, oo)) of g (x) and / (x), respectively. Also, assume that both the harmonic 
extensions vanish for \x\ ,\y\ — > oo, then 



/ V n+ iu(x,y) ■V n+ iv(x,y)dxdy = / V n+r (u (x, y) V n+ iv (x, y)) dxdy 
Jc Jc 

u(x,y)V n+ iv (x,y) dxdy 

da(121) 



dC 

oy 



Here, dC indicates the border of C. To obtain this identity we have written 
the integral of the (n + l)-dimensional divergence (|12ip as a surface integral 
over the component of the vector field normal to the surface dC which is 
taken as a finite (n + l)-dimensional cube with a face at y = 0. Then we 
have made the 2n + 1 faces other than y = to go to infinity, where functions 
u and v vanish. 

Hence, for the harmonic extensions u,v to C = M n x (0, oo) of g,f, re- 
spectively, we have 



,(„rtvs tl .(„»)**-^(„rtvs +1 . fe ,)** 

Jc 
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or, using f|121j) . 

dx = 



d d 
u [x, y) —v (x, y) -v (x, y) —u (x, y) 

In terms of g (x) and / (x) this identity reads 



y=0 



,<.)^-*&/m)*-o 

or, equivalently 

g{x)Tf{x)dx = f(x)Tg(x)dx (122) 



Concerning the case in which fermions are coupled to the U(l) gauge 
field, we introduce the operator T[A] defined as 

T[A] = (-A[^]A[^]) 1/2 (123) 

with the covariant derivative given by 

Di = idi + eAi (124) 

In the case of d = 1 space dimensions this simplifies to 

T[A] = (-D 1 [A]DM]) 1/2 (125) 

The covariant derivative D\ [A] can be written as 

Z?i[A]=exp^e J Axdy^j (^i-^J exp (-ie J A x dy\ (126) 

which in terms of a Wilson line 

W[Ax] = exp (^-ie J Axdy^j (127) 



takes the compact form 



D X [A] = W^A^i^-W^] (128) 



One then has 



{DM]? = W- 1 ^} (i-fLYwlAi] = W'M^WiAr] (129) 



or 



(D^A}) 2 = W-^AtlTWiA^W^lA^TWiAt] (130) 

so that 

T[Ai] = W-^A^TWiAr] (131) 
showing that T[Ai] is Hermitian. 
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